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We investigate the effects of quenched disorder on the coherent response in a driven system of 
coupled oscillators. In particular, the interplay between quenched noise and periodic driving is 
explored, with particular attention to the possibility of resonance. The phase velocity is examined 
as the response of the system; revealed is the enhancement of the fraction of oscillators locked to 
the periodic driving, displaying resonance behavior. It is thus concluded that resonance behavior 
may also be induced by quenched disorder which does not have time-dependent correlations. 
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In recent years, the stochastic resonance (SR) has 
drawn much attention, which stands for the phenomena 
that the response of a system to a periodic driving force is 
enhanced by an appropriate amount of noise rather than 
suppressed 1, 2, 3]. Those SR phenomena observed in 
various real systems 5^ _6] are known to occur through 
the cooperative interplay between the noise and the ex- 
ternal driving force, and require three ingredients: (i) 
an energetic activation barrier; (ii) a weak input such 
as periodic driving; (iii) a random noise 3]. In par- 
ticular, the SR may be understood in terms of match- 
ing of two time scales, the period of the driving force 
and the inverse of the Kramers hopping rate associated 
with the random noise. Note that the noise here usu- 
ally has time-dependent correlations. Namely, the cor- 
relation between two noise forces rji and rjj is given by 
{r]i{t)rjj [t')) = 2T6ij6{t—t'), where T represents the noise 
strength. Here naturally arises a question how does the 
SR behavior depends on such characteristics of the noise. 
In particular, one may ask whether such resonance be- 
havior appears even for the quenched noise without time- 
dependent correlations. 

To resolve this, we consider the system of coupled os- 
cillators described by 



= uJi - —^^sin{ipi ~ ipj) + licosilt, (1) 
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where (pi represents the phase of the ith limit-cycle os- 
cillator. The first term un on the right-hand side de- 
notes the intrinsic frequency of the zth oscillator, as- 
sumed to be distributed according to the Gaussian dis- 
tribution (/(w) with zero mean ((w) = 0) and correlations 



D5ij. Note that correlations between w^'s are 



not time-dependent but quenched in time, in sharp con- 
trast with the usual thermal noise. The second term on 
the right-hand side describes the coupling between os- 
cillators, whereas the third term represents the external 
periodic driving with the driving strength Ii chosen from 
a certain distribution function /(/). In the absence of 
the external driving [Ii — 0), Eq. (Q) describes the well- 
known Kuramoto model, where the synchronization phe- 



nomena have been investigated extensively . The scat- 
tered intrinsic frequency or quenched noise competes 
with the coupling strength K in the system: When the 
coupling is strong enough to overcome the scatteredness 
of intrinsic frequencies, synchronization emerges. Such 
synchronization behavior may be explored by measuring 
the order parameter 



(2) 



where nonzero A implies the emergence of synchroniza- 
tion. It is known that the synchronization behavior of the 
system is characterized by A [K — with /3 = 1/2 
near the critical coupling strength Kc = 2/'Kg{0) 7]. The 
Kuramoto model has been also extended by means of in- 
troducing a constant in the argument of sine coupling |^ . 
Such extended model is also studied in Ref. '9|, where 
frequency locking in the system of Josephson arrays is 
investigated. 

Meanwhile, when the external periodic driving comes 
into system, the synchronization behavior has been ob- 
served to appear periodically 0|. However, the possi- 
bility of the resonance phenomena, which may occur via 
the cooperative interplay between the external driving 
and the quenched noise, has not been investigated. In 
this work, we examine the linear response of the system, 
with particular attention to the possibility of the reso- 
nance behavior due to the quenched disorder. 

We first investigate analytically the dynamics of the 
system governed by Eq. The order parameter A 

defined in Eq. Q allows us to reduce Eq. ^ into a single 
equation 



(p = uj — KA sin (f + I cos fit, 



(3) 



where A is to be determined by imposing self-consistency, 
and indices are suppressed for simplicity. Equation Q 
reminds one of the single resistively shunted Josephson 
junction under combined direct and alternating currents. 
It is well known that such a system can be locked to the 
external driving, which is characterized by the Shapiro 
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with n integer, where ip (= t;) is the time-averaged phase 
velocity. Such mode locking features suggest the ansatz 

oo 

ip = Lpn + nilt + Ap sm{pnt + Up) (5) 
p=i 

for the locked phase of the oscillator on the nth step. 
The external periodic driving in the system also leads us 
to expect periodic synchronization [lOj, where the order 
parameter A is decomposed as 



A = Ao + ^ A, cos{qnt + /3,). 
Inserting Eqs. and © into Eq. ||2Jl, we obtain 

oo 

nil + pilAq cos{prtt + ap) — I cos Q,t 

oc 

Ao + ^Aq cosiqflt + Pq) 

9=1 

/ OO OC \ 
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= uj-K 
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where Ji^{x) is the ^^th Bessel function and $ = (^o + 
Yl,T=i ^sOis + + sis)^t. The integers £s satisfying 
sts = ~" contribute to the dc component in <I>: 

nVl^oj- KaJyYY^ JiMr)] sin$ + C'(i^Ai), (8) 

\r=l ) 

where $ = (^sq + X^^i ^s^s ^-nd the prime in the summa- 
tion stands for the constraint sis — ~n. This gives 
an estimation of the dc driving strength lo corresponding 
to the integer locking. The amplitude Ap and phase ap of 
the ac component with frequency pVl can be determined 
from the equation 

— pflAp cos{pQt + ap) + Sp il cos pflt 

= KAa lf[^Je+ (Ar)] sin l^o+P^t + f^ ^+a, ) 



+ KAo Y[ J2 -^e- (Ar) sin iipo-pnt + Y, C". 

\r=l g- / \ s=l ) 

+ KAp cos{pnt + /3p) ( n X! ^^"-n * 

with integers and f ~ satisfying ^^t = p — n 

and '5^7 ~ ~P ~ respectively. When KAq is 



sufficiently small compared with the driving amplitude, 
Eq. ^ with p = 1 yields Ai and ai to the zeroth order in 
KAq: Ai = I/ft and ai = 0. It is observed numerically 
that the dc component Aq of the order parameter is much 
larger than higher-order (ac) components [Aq ^ A^ > 
1)], which allows the expansion of the locked phase to 
the zeroth order in KAq. This gives the (locked) phase 
of the oscillator on the nth step as 

ip = Lpo+ nnt + ^ sin Qt + 0{KAoI/2n'^), (10) 

which yields the range of locked oscillators: 

uj-nQ^ if Ao(-l)" J„(//17) sin 00 + O(A'Ai). (11) 

It implies that the oscillators in the range 

nfl - KAo\JniI/^)\ <LU<nn + isTAol J„(//f7)| (12) 

display the locking behavior (v /f2 — n) , with the higher- 
order terms of 0{K Ai) neglected. 

The quenched disorder uj and the driving amplitude / 
are chosen from the distribution function g{Lu) and /(/), 
respectively, which yields the fraction r„ of the oscillators 
locked to the nth step: 



f{l)dl / g{uj)duj. 



(13) 



For the Gaussian distribution g{uj) = (l/\/27rD)e"" /^-^ 
and the delta function one /(/) = (1/2) [(5(7 - /q) + 5(7 + 
Iq)], the fraction r„ in Eq. (|13|l is given by 
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where erf(a:) denotes the error function. To estimate the 
behavior of the fraction r„ as the variance D varies, we 
should know the behavior of the dc component Aq as a 
function of the variance D. Note that the self-consistency 
equation for the order parameter gives the behavior of Ag 
only near the critical point. To obtain Aq in the whole 
range of the disorder strength, we resort to numerical 
simulations and integrate Eq. (Q) via Heun's method (l3| 
with the discrete time step St = 0.01. While the equa- 
tions of motion are integrated for Nt = 8 x time 
steps, the data from the first Nt/2 steps are discarded in 
measuring quantities of interest. The system size A'" has 
been considered up to iV = 20000, so that no appreciable 
size-dependence is observed. The driving amplitude Iq in 
the distribution /(/) and the driving frequency Q have 
been chosen to be Iq = 0.8 and O = 1.0, 1.2, and 1.4, 
respectively; the coupling strength K has been set equal 
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to unity for convenience. We measure the order parame- 
ter A, and obtain the component Aq by taking the time 
average. 

The do component Aq is shown in the inset of Fig. Q] 
as a function of the disorder strength £>, displaying the 
monotonicaUy decreasing behavior. Using this, we com- 
pute the fraction r±i of the oscihators which are locked 
to the first step since r±i is most dominant over other 
components [r^(f > 2)]. Figure displays the total frac- 
tion ri + r_i for K = 1.0, Iq = 0.8, and O = 1.0, 1.2, 1.4 
versus the disorder strength D. It is found that the frac- 
tion first increases with D, which implies that a larger 
number of oscillators tends to follow the driving force as 
the quenched disorder becomes stronger. Remarkably as 
the disorder is increased further, the fraction reaches the 
maximum and begins to decrease. For example, the op- 
timal disorder strength is observed to be Dm ~ 0.17 for 
il — 1.0. Such behavior of the fraction r±i is reminiscent 
of the SR, which is known to occur through the coop- 
erative interplay between the external periodic driving 
force and the random noise. Note that in sharp contrast 
to the random noise in the conventional SR, the noise 
in the system governed by Eq. is quenched with no 
time-dependent correlations. Namely, here the quenched 
disorder enhances the coherent response of the system. 
It is also shown that as the driving frequency fl is raised, 
the fraction r±i of locked oscillators diminishes while the 
optimal strength Dm shifts to larger values. These ten- 
dencies reflect that the oscillators are reluctant to follow 
the driving which changes too fast, and that the disor- 
der strength D should be enlarged to fit the high driving 
frequency. Another point is that the fraction vq of the 
oscillators locked to the n = step does not exhibit such 
resonance behavior; it rather displays the monotonic de- 
creasing behavior, which is quite similar to that of the 
dc component Aq of the phase order parameter. Higher 
fractions for n > 2 also show the resonance behavior, al- 
though the magnitude is too small to be clearly discrim- 
inated. We have also investigated the role of coupling on 
the enhancement of coherent response. As the coupling 
strength K decreases, such resonance behavior is found 
to be suppressed. 

To confirm such resonance behavior, we now perform 
numerical simulations. We obtain the probability distri- 
bution P{v) for the time-averaged phase velocity v un- 
der the same conditions as before, varying the disorder 
strength D. The data are then averaged over 100 inde- 
pendent sets of {uJi}. The inset of Fig. El displays the 
probability distribution for f2 = 1.0 and D — 0.20, where 
three sharp peaks appear at w = (not fully shown) and 
±0. The two peaks at the driving frequency (±f2) man- 
ifest that some fraction of the oscillators follow the ex- 
ternal driving force, locked to the latter. Note that such 
peaks may also appear at higher frequencies nil (|n| > 2), 
although they are too small to be manifestly shown. 

In the absence of the quenched disorder {D — 0), those 
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FIG. 1: Total fraction ri + r_i of the oscillators locked to 
the first step (n = ±1) versus the disorder strength D for 
K = 1.0, Jo = 0.8, and fl = 1.0, 1.2, and 1.4. Inset: dc 
component Aq in the system of size = 20000 at various 
driving frequencies. 



peaks a.t V = ibfl do not emerge, and only the delta peak 
appears at u = 0. As the disorder strength is increased 
further, on the other hand, peaks show up at the driv- 
ing frequency, with the height growing. The peak height 
h±i at the driving frequency may be regarded as a good 
indicator which describes the response of the system to 
the external driving. 

We numerically obtain the peak height, subtracting 
the background noise given by the average of ten data 
points around the peak. Figure|21displays the total height 
hi + h^i versus the disorder strength D in the system 
of size N = 20000, with the driving frequency varied 
from 17 = 1.0 to = 1.4. We have considered the sys- 
tem size up to = 40000, where no appreciable size- 
dependence is observed. The total height is found to 
first increase with the disorder strength D and reaches 
its maximum at a finite value of the disorder, displaying 
quite a similar feature to the fraction r±i [see Fig. The 
enhancement of the height indicates that the quenched 
disorder actually increases the number of the oscillators 
locked to the external driving. To our knowledge, such 
enhancement induced by the quenched disorder without 
time-dependent correlations has not been addressed be- 
fore. It is observed in Fig. El that the increase of the 
driving frequency tends to suppress the total height, and 
to shift the optimal disorder strength to larger values, 
which may be related with the intrinsic time scale of the 
system. Note that the conventional SR phenomena have 
been known to occur when the Kramers hopping rate and 
the external driving frequency match each other. To see 
such time scale matching in the quenched-noise-induced 
resonance, we investigate the relaxation dynamics of the 
system, and probe the time evolution of the renormalized 
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FIG. 2: Total height hi + h-i of the peak at u = ±57 in the 
probability distribution versus the disorder strength D. In- 
set: probability distribution P{v) of the time-averaged phase 
velocity v ior K = 1.0, Jo = 0.8, Q. = 1.0, and D = 0.20. 
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FIG. 3: Inverse relaxation time vs the disorder strength 
D. The crossing points with the lines 1.0/27r, 1.2/27r, and 
1.4/27r yield the values of Dm at the corresponding driving 
frequency. 



synchronization order parameter 

A(t) - Ac 



A{t) 



A(0) - Ae 



(15) 



where Acq and A(0) represent the equiUbrium value and 
the initial one, respectively. The renormahzed order pa- 
rameter A(t) is thus expected to decay from unity (t = 0) 
to zero (t — > oo): A{t) exp(— f/r). 

Figure 13 shows the inverse relaxation time r^^ versus 
the disorder strength D for Iq — 0.8 and fl — 1.0, 1.2, 
and 1.4. The intrinsic time scale r of the system is ob- 
served to increase indefinitely as the disorder strength 
D approaches the critical value Dc beyond which syn- 
chronization disappears. The horizontal lines in Fig. El 
describe the time-scale matching conditions at various 



driving frequencies: the value of D at the crossing point 
corresponds to Dm at which the enhancement of the re- 
sponse is maximized. For example, for ft — 1.0, 1.2, and 
1.4, the crossing points lead to Dm ~ 0.19, 0.21, and 
0.22, respectively, which are consistent with the values 
obtained both analytically and numerically [see Figs. 
andU]. We thus conclude that quenched disorder may 
also enhance the coherent response of the system via the 
mechanism of the time-scale matching. 

In summary, we have explored the effects of quenched 
disorder on the coherent response in the driven system 
of coupled oscillators. We have investigated the interplay 
between quenched disorder and external periodic driving, 
with particular attention to the possibility of resonance 
behavior. The phase velocity is probed as the response of 
the system; revealed is the enhancement of the fraction of 
the oscillators locked to the external driving, exhibiting 
resonance behavior. This provides the first observation 
of the resonance behavior induced by quenched disorder. 
In a biological system such as synchronous fireflies, dif- 
ferent firing frequencies of fireflies may be regarded as the 
quenched disorder. Our results are applicable to expect 
those different firing frequencies (instead of same ones) 
may enhance the coherent response of the system. 

We thank M.Y. Choi, H. Park, and B.J. Kim for useful 
discussions. 
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